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ABSTRACT 


A heat-transfer analysis was made for a coaxial-flow nuclear rocket. The analysis 
considers one -dimensional radial transfer of heat, by both radiation and turbulent mix- 
ing from a centrally located fissioning gas to a coaxially flowing propellant. The results 
compare the effective thermal conductivities for radiation and turbulent mixing. Also 
included are typical radial-temperature profiles at various axial locations and average 
fuel temperatures showing the effect of neglecting turbulent mixing. Turbulent mixing 
had a large effect on local temperature in regions where the fuel density was smallest. 
Therefore, the effect of turbulent mixing on the average fuel temperature was small. 



EFFECT OF TURBULENT MIXING ON AVERAGE FUEL TEMPERATURES 
IN A GAS-CORE NUCLEAR ROCKET ENGINE 
by Albert F. Kascak and Annie J. Easley 
Lewis Research Center 

SUMMARY 

Heat transfer in a coaxial-flow, gas-fueled nuclear rocket with a specific impulse of 
1500 seconds and a thrust between 0. 445 and 4. 45 meganewtons (0. 1 and 1. 0 Mlb) was 
analyzed. In this type of rocket, a low-velocity, fissioning, gaseous fuel is injected so 
as to flow axially down the center of a cylindrical reactor cavity. A high-velocity, hy- 
drogen propellant flows coaxially around the fuel. Heat is generated by the fissioning 
fuel and is transported to the propellant. Past analyses have assumed radiation to be the 
principal mechanism of heat transfer. Other analyses have shown turbulent mixing to be 
the dominant means of mass transfer. The object of this study was to determine the ef- 
fect of turbulent mixing on the transport of heat. 

The present analysis considers one -dimensional radial transfer of heat by both radi- 
ation and turbulent mixing. Effective thermal conductivities for radiation and for turbu- 
lent mixing are used. The effective radiation conductivity is based on the Rosseland dif- 
fusion formulation of radiative transport, and the turbulent conductivity is based on an 
empirical eddy viscosity model. 

The results of this study include the following: 

1. Eddy viscosities obtained by scaling isothermal fluid dynamic results 

2. Comparison of the effective conductivities for radiation and turbulent mixing 

3. Typical radial-temperature profiles at selected axial locations 

4. Average fuel temperatures showing the effect of turbulent mixing 

These results showed that the change in the average fuel temperature was less than 
1 percent when the heat transferred by the turbulent mixing was included. Turbulent 
mixing had a large effect on temperatures in regions where the fuel density was smallest. 
This conclusion is valid unless eddy viscosities have been underestimated by at least 
several orders of magnitude. Although these conclusions were drawn for a particular- 
reactor geometry and flow distribution, changes in these factors would not cause order- 
of-magnitude changes in the eddy viscosities. 



INTRODUCTION 


The gas-fueled nuclear rocket is a proposed propulsion concept that has high specific 
impulse and relatively high thrust. One configuration, illustrated in figure 1(a), is known 
as the coaxial -flow,' gas-fueled nuclear rocket. In this concept, a low -velocity, fission- 
ing gas fuel flows along the axis of a cylindrical cavity. The propellant flows coaxially at 
high velocity around the central fuel region. The objectives of using this concept are to 
transfer energy from the hot fuel to the colder propellant and to minimize erosion of the 
fuel by the propellant. 

To minimize fuel-propellant mixing, an intermediate -velocity buffer layer is placed 
adjacent to the fuel. The buffer layer may be either hydrogen or some other gas whose 
molecular weight is between that of the fuel and the propellant. The propellant to be 
heated is seeded with opaque solid particles to increase its opacity; it flows between the 
fuel-buffer layer and the porous wall. Unseeded propellant provides transpiration cool- 
ing of the porous wall. Reference 1 recommends bypassing part of the propellant to in- 
crease the fueled region and to reduce the seeded propellant velocity and thereby reduce 
fuel-propellant mixing. Due to mixing, these flow regions are not distinguishable at all 
points in the reactor cavity. 

Past analyses have assumed radiation to be the dominant mode of heat transfer be- 
cause of the high fuel and propellant temperatures. But reference 2 has shown that tur- 
bulent mixing is the dominant means of mass transfer in the fuel -propellant mixing pro- 
cess. Heat transfer by turbulent mixing could be important where temperature gradients 
or shear stresses, or both, are large. Three such regions exist in the coaxial -flow 
rocket. The first is near the inlet where large temperature gradients exist between the 
fuel and the propellant. A second is along the walls of the cavity and the nozzle, and a 
third is in the fuel-propellant jet mixing region. In the last two regions shear stresses 
would be large. 

Three questions are to be answered in this report: 

(1) Does turbulent mixing transfer a significant fraction of the fission power from 
the fuel to the surrounding propellant ? 

(2) Is the effect of turbulent mixing dependent on engine thrust or reactor pressure? 

(3) Is this estimate of the effect of turbulent mixing or heat transfer applicable to a 
wide range of operating conditions? 
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(b) Representative cross section showing pertinent dimensions. (At z = 0, r o (0) = n. ) 

y J 

Figure 1. - Coaxial-flow, gas-fueled nuclear rocket. 



ROCKET OPERATING CONDITIONS 


In order to perform the analysis, a suitable set of operating conditions must be se- 
lected. Reference 3 lists the following representative values for a first-generation gas- 
core nuclear rocket: 


Specific impulse, I , sec . . . . 
sp 

Fuel- to total-flow -rate ratio, r ^ 
Thrust, F, MN (Mlb) 


1500 

1/31 

0.445 to 4.45 (0. 1 to 1.0) 


From these values, the following mass flow rates may be calculated: 


Exhaust, m^ = F/(gI gp ), kg/sec 34 to 340 

Fuel, m f = r? f m ex , kg/sec 1.1 to 11.0 

Propellant, m p = (1 - Vf)™ ex , kg/sec 32.9 to 329 


When a nozzle efficiency r i _ of 90 percent is assumed, the exhaust enthalpy is 

noz 


h 


ex 



1.33X10 8 J/kg 


( 1 ) 


Reference 4 studied the use of a buffer layer to reduce fuel-propellant mixing and 
concluded that mixing was minimized for Up/U^ = 5 and approximately equal flow areas 
^r^ = r. + 1/2 (r w - r^. For this study two-thirds of the propellant was bypassed, one- 
sixth of the propellant was used for wall cooling, and one -sixth of the propellant was in- 
jected coaxially to the fuel as the buffer layer and the seeded propellant. 

Reference 3 also indicated that both the diameter and length of the reactor cavity 
should be approximately 3 meters. A cavity diameter of 3 meters and a cavity length of 
2. 1 meters were chosen for this analysis. The diameter of the fuel region at the inlet 
was assumed to be 2. 1 meters. Therefore, the maximum fuel volume was 50 percent of 
the total cavity volume with the intermediate bypass mixing region neglected. References 
5 to 7 indicate that the critical mass is between 10 and 100 kilograms. 

When only radiation is considered, reference 8 found the average fuel temperature 

to be about 50 000 K. Therefore, the fuel is ionized about four times. According to the 

7 8 

perfect gas law, the operating pressure is between 10 and 10 newtons per square 

o o 

meter (10 and 10 atm). (If not specified differently, the operating pressure was as- 
sumed to be 5. 05x10 N/m (500 atm). ) The average propellant temperature can be cal- 
culated by assuming that all the heat is transferred from the fuel to the propellant be- 
tween the fueled region and the porous wall. Because this propellant was assumed to be 
one -third of the total amount of propellant, the enthalpy at the exit of the cavity region 
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(prior to mixing with the bypassed coolant) must be three times the exhaust enthalpy. 

For a rocket with a specific impulse of 1500 seconds using hydrogen as the propellant, 
this exhaust enthalpy corresponds to a propellant temperature at the exit of the cavity re- 
gion of about 10 000 K. The average propellant temperature in the cavity is about one- 
half this temperature, or 5000 K. When this temperature and the average fuel tempera- 
ture are used, the average fuel density is three times the average propellant density. 

The Reynolds number of the injected fuel is defined as 


p.U.r. m f 
Re. = J . 3 I = — L 


3 




3 


Ir i u j 


or, in terms of thrust, 



7] f F 


( 2 ) 


( 3 ) 


For the range of thrusts considered, the Reynolds number for the fuel jet is between 
7.67X10 3 and 7.67X10 4 . 

The energy generated in the fuel and transferred directly to the propellant is between 
3. 97xl0 3 and 3. 97xl0 4 megawatts where 

Qt = “p h ex “ Qmod 

The energy dissipated in the moderator is taken to be 0. 10 Qj. as given in refer- 

ence 3. 

The authors believe these operating conditions to be representative of first- 
generation gas-core nuclear rockets. Characteristics of such a rocket should be com- 
mon to most gas-core nuclear rockets. Therefore, the conclusions drawn should be in- 
sensitive to the particular numerical choices made. A summary of the important oper- 
ating conditions is given in table I. 


ANALYSIS 

The analysis is divided into three parts. The first section discusses the formulation 
of approximate laws for determining the properties of mixtures from the properties of 
their constituents. Second, a fluid dynamics analysis is used to determine the amount of 
turbulent mixing. Third, formulations of the heat -transfer equation are presented. 
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Mixture Properties 


Approximate mixing laws must be known so that the properties of the fuel-propellant 
mixture may be calculated from known properties of the fuel and propellant. The proper- 
ties needed for this analysis are the Rosseland mean absorption coefficient, the specific 
heat, the viscosity, and the fuel density. 

The Rosseland mean absorption coefficient is based on the average thermal radiation 
cross section of a particle; thus, its mixing law is related to the atom fraction, or 

a R = a R, f x + a R, p^ 1 " x ) 

Since specific heat is defined on a mass basis, its mixing law is related to the mass 
fraction so that 


c p = c p,f y * c p,p ( 1 _y) (6) 

Viscosity is a transport property and therefore is not the linear sum of the viscosity of 
the constituents. Reference 9 gives the following mixing law for viscosity. 


where 


and 


ju = 


xfi f 


(1 - x)p 


P 


x + (1 -x)<p fjP X<p p>f +(l-x) 



By the perfect gas law, the fuel density is 


p f = 


M f P f 

RT 


(7) 


( 8 ) 


(9) 


( 10 ) 


6 



and by the definition of the mole fraction, the fuel partial pressure is 


P f = xP t 


( 11 ) 


the mass and mole fractions and the average molecular weight are defined as follows: 


y = 


Pf 

p f + p p 


( 12 ) 


X = 


p f + p p 


y 


/ Mp\ 

y + ( — ) (i - y) 

K/ 


(13) 


— . 
M = 1 


0! 


S.n 


1 CH, 1 


(14) 


where n . is the number density of the pure a substance of the i 1 n specie. These 

Cz, 1 

mixing laws are approximate, and they do not allow for a change in individual specie con- 
centration when the individual substances are mixed. The authors assume these approx- 
imations to be sufficiently accurate for this analysis. 

The Rosseland mean absorption coefficient for the propellant is given in refer- 
ence 10. The data from which the specific heat and molecular weight were calculated 
are also given in this reference. The viscosity of the propellant was taken from refer- 
ences 11 and 12, and its molecular conductivity was obtained from reference 13. 

The viscosity of the fuel from reference 11 was extrapolated to lower temperatures. 
The Rosseland mean absorption coefficient of the fuel is given in reference 14. The spe- 
cific heat and molecular weight of the fuel were computed by the method given in refer- 
ence 15 after modification to include charged particles. Fuel ionization potentials used 
in this calculation were provided by J. T. Waber, D. Lieberman, and D. T. Cromer of 
Los Alamos Scientific Laboratory. These potentials are presented in reference 8. 
References 16 and 17 describe related work. 


Fluid Dynamics 

In order to analyze the effect of turbulent mixing on heat transfer, the magnitude of 
the mixing must be determined from a fluid dynamics analysis. Reference 18 describes 
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a computer code for determining the isothermal, turbulent-jet mixing of two gases. This 
code was used to compute the axial mass flow rate m, the eddy viscosity ep, and the 
containment factor I. Since the code considers only isothermal mixing, its results must 
be scaled by forming dimensionless groups which the authors assume will not be temper- 
ature dependent. This scaling process is only approximate since the energy equation was 
not solved by the code. 

Two empirical eddy diffusivity correlations based on a characteristic length are used 
in the code. Before the fuel jet centerline velocity begins to change, a wake flow corre- 
lation is used in which the axial distance from the inlet is the characteristic length. Fur- 
ther downstream where the mixing is axisymmetric, a fully developed correlation is used. 
This correlation is based on the half radius, which is the point at which the velocity 
equals the average of the fuel centerline and propellant velocities. 

The containment factor is defined as 


I(z) =— L- 

„ „3 



pyr’ dr’ dz’ 


(15) 


where I(z) is the amount of fuel contained in the cavity to axial position z. 

The velocity ratios U^/Uj = 5 and Up/Uj = 25 and the density ratio Pp/Pj = 3, 
which were evaluated in the section ROCKET OPERATING CONDITIONS, were used as 
part of the code input. The mass fraction y, the dimensionless mass flow rate pU/PjUj, 
and the dimensionless turbulent viscosity e P/ r jUjPj were taken from the code output. 

The dimensionless parameters y, pU/p.U-, and (ep/p)/(r.U-P-/p.) = e + /Re. were as- 

J J J J J J J 

sumed by the authors to be independent of temperature. 


Heat Transfer 

The relations governing radiant heat transfer in the coaxial-flow, gas-core nuclear 
rocket were developed in reference 8. This work must be extended to include turbulent 
mixing. 

The analysis considers one -dimensional radial transfer of heat due to both thermal 
radiation and turbulent mixing. The formula 

q=-K^I (16) 

3r 

allows the definition of an effective thermal conductivity K where 
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K = K rad + ^urb 


(17) 


The effective radiation conductivity K ra( j is based on the Rosseland diffusion formu- 
lation of radiative transport. The effective turbulent conductivity is based on an 

empirically defined eddy diffusivity. 

The effective radiation conductivity is inversely proportional to the Rosseland mean 
absorption coefficient and directly proportional to the temperature cubed. The resulting 
equation is 


K. 


rad 


16 ffT‘ 


3 a 


(18) 


R 


The effective turbulent conductivity is equal to the eddy viscosity ep times the specific 
heat so that 


Kturb = 

where 


(19) 


e 


+ 


£ 

V 


€ 




m 


The local heat-source strength is proportional to the local fuel density when the 
neutron flux is assumed to be uniform with respect to position and energy. When 
Q 1 "/ Pj is the coefficient of proportionality, the equation for local heat source strength 
is 


Q 


m t 



( 20 ) 


Performing a heat balance on the interior regions of the core results in a relation be- 
tween the heat-source strength and the heat flux. 


27rrq(r) = 2 ir 


f 

• 4 ) 


r , Q ,M (r , )dr’ 


( 21 ) 


The average fuel density is defined as the density of pure fuel at the system total 
pressure and at an average temperature defined implicitly by 
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mil in win mi in iniii i mm i Him 


in i 


/*P r 

in 6 P f ( T, y)r* dr' = 2v e p f (T, y)r* dr' 


( 22 ) 


where y is the mass fraction at any radius and the fuel-propellant boundary r g is the 
radius at which y = 0. 01. The average fuel density is then pif = P^(t, 1. o). The average 
power level Q M * is related to heat flux at the fuel -propellant boundary and the average 
temperature. This relation can be found by combining equations (20) to (22) 


Q’” = 


■eVf 


P f (T, y)r* dr* 


(23) 


The blackbody radiating temperature at the fuel-propellant boundary is defined by 


q 


e 



(24) 


The actual temperature is some multiple of the boundary blackbody radiating temperature 

T e = tjT* (25) 

where t] = 1 (based on ref. 8). The radial heat flux is given by 


q !<1_ 

6 2irr e r j I ( L ) d( z Aj) 


(26) 


Equations (16) to (25) are a set of equations which describe the heat -transfer pro- 
cess. Appendix A describes an iterative Runge-Kutta technique for solving these equa- 
tions. 


RESULTS AND DISCUSSION 

The results obtained are discussed in this section. First, the gas mixing patterns 
and the eddy viscosities used are discussed. Next, the effective thermal conductivities 
due to both turbulence and thermal radiation are presented and compared. Finally, fuel 
temperatures in an engine are calculated first with turbulence neglected, and then with it 
included. 

Figure 2(a) shows the dimensionless axial mass flow rate PU/PjUj as a function of 
the dimensionless radius r/r. for various values of the dimensionless length z/r,. The 

J J 
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Mass fraction of fuel, 


! 



0 .4 .8 1.2 1.6 


Nondimensional radius, r/r j 

(a) Nondimensional mass flow rate. High- 
velocity outer stream first accelerates 
buffer layer then accelerates fueled region. 



0 .4 .8 1.2 1.6 0 .4 .8 1.2 1.6 


Nondimensional radius, r/r j 

(b) Mass fraction of fuel. Fuel first diffuses (c) Normalized eddy viscosity. Eddy viscosity 

into buffer layer then is eroded away by at inlet is low and then increases with 

high-velocity propellant. increasing axial position. Eddy viscosity 

is lower in fuel than in propellant. 

Figure 2. - Nondimensional mass flow rate, mass fraction of fuel, and normalized eddy viscosity as 
functions of nondimensional radius for various axial positions. 
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fuel velocity was assumed to change instantly at the point z = 0, where the fuel enters 
the fuel region. At z/r^ = 0. 5 the mass flow rate of the buffer layer is about one-half 
that of the free stream. The mass flow rate of the buffer layer approaches that of the 
free stream for z/r^ = 1, and, for z/r^ > 1, the fuel is accelerated by the free stream. 

Figure 2(b) shows the fuel mass fraction as a function of r/r^ for various values of 
z/r j. When a step change in the fuel mass fraction at z = 0 is assumed, the fuel radius 
is a maximum at z/r^. = 0. 5. For z/r^ > 1 the slower moving buffer layer no longer 
exists. Therefore, the fuel diffuses into the high-velocity propellant and is swept away. 

Figure 2(c) shows the normalized eddy viscosity e + /Rej as a function of r/r^ for 
various z / v y At the inlet e + /rtej is assumed to be zero. In general, e + /Rej is almost 
constant in the propellant region and decays in the fuel region. The magnitude of e + /Rej 
increases with increasing z / r y 



Figure 3. - Containment factor as function of 
nondimensional axial position. Near inlet 
there is very little mixing. Mixing process 
increases downstream, causing containment 
to fall below perfect containment line. 


Figure 3 shows the containment factor I as a function of z / r y E there were per- 
fect containment (no mixing), the slope of the line would be 1/2. Since there is mixing, 
the actual line falls below the perfect containment line. The slope of this line is less 
than 1/2 and decreases with increasing z/r.. Physically, the mixing rate starts at low 
values and increases with increasing z/r.. 

Figures 2 and 3 show that the majority of the mixing occurs in the downstream half 
of the reactor. For a thrust of 2. 5 meganewtons (0. 5 Mlb), the Reynolds number for the 
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fuel jet is 38 300. For this Reynolds number, figure 2(c) shows that the dimensionless 
eddy viscosity e + is between 200 and 2000 for z/r^ > 1. For the constant dimension- 
less eddy viscosity calculations, a value of 1000 was used. 

Figure 4(a) shows the effective turbulent mixing and radiation conductivities for both 
the fuel and the propellant and the molecular conductivity for the propellant. The molec- 
ular conductivity of the fuel is expected to be one to two orders of magnitude less than 
that of the propellant because, although the electron densities are comparable, the spe- 
cific heats differ by one or two orders of magnitude. The molecular conductivities are 



(a) Conductivity for pure constituents for 
various modes of heat transfer. Pro- 
pellant conductivities are several 
orders of magnitude higher than re- 
spective fuel conductivities. At high 
temperatures radiation conductivities 
are several orders of magnitude higher 
than turbulent conductivities. 


(b) Combined conductivities for various 
mass fractions. Turbulence is im- 
portant below temperatures of about 
30 000 K. In general, mixtures of 
fuel and propellant tend to have con- 
ductivities similar to propellant. 


Figure 4. - Conductivity as function of temperature at pressure of 5.05x10? newtons per 
square meter (500 atm). 
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composed mainly of an electron conductivity and a reaction conductivity which is propor- 
tional to the specific heat. In this temperature range, the reaction conductivity is domi- 
nant. The effective radiation and turbulent mixing conductivities for the propellant are 
also one to two orders of magnitude higher than the corresponding fuel conductivities. 
Radiation is the dominant mode of heat transfer for temperatures greater than 25 000 K, 
and for temperatures greater than 50 000 K radiation is several orders of magnitude 
greater than turbulent mixing. Below 20 000 K turbulent mixing is the dominant mode of 
heat transfer. In the range of temperatures shown, molecular conduction is negligible. 

By using the mixing laws given in the section ANALYSIS, the total effective conduc- 
tivities were calculated as a function of temperature for various mass fractions. These 
effective conductivities are shown in figure 4(b) with and without turbulent mixing. At 
temperatures below 30 000 K turbulent mixing greatly increases the total effective con- 
ductivity, while above 30 000 K the effect is small. Figures 4(a) and (b) also show that 
the heat -transfer process is a variable -properties problem, so that the total conductivity 
is a strong function of the propellant properties. 

For the operating condition given herein, an average radial temperature gradient of 
42 500 K per meter was calculated. When a reactor power of 20 000 megawatts (2. 5-MN 
thrust) and the heat -transfer area of a cylinder whose radius and length are 2. 1 meters 
are used, an average heat flux of 1. 42x10 megawatts per square meter was calculated. 
The average total conductivity must therefore be about 33 kilowatts per meter per Kelvin. 
This conductivity corresponds to a fuel temperature of about 63 000 K or a propellant 
temperature of 35 000 K. Figure 4(b) shows that at these temperatures turbulent mixing 
has a negligible effect on the total conductivity. Figure 4(a) shows that for turbulent mix- 
ing to be dominant at these temperatures, the estimate of the eddy viscosity must be in 
error by several orders of magnitude. 

Figures 5(a) and (b) show the effect of turbulent mixing on the radial temperature 

7 

profiles at various axial locations for a rocket operating at a pressure of 5. 05x10 new- 
tons per square meter (500 atm), a specific impulse of 1500 seconds, and thrusts of 
0.445 and 4.45 meganewtons (0. 1 and 1.0 Mlb), respectively. Turbulent mixing reduces 
the temperature in the downstream jet mixing region by about 30 percent for a thrust of 
0. 445 meganewton (0. 1 Mlb). At a thrust of 4. 45 meganewtons (1. 0 Mlb) the change in 
temperature was about 90 percent. The effect near the inlet or centerline is less than 
1 percent. Comparison of figures 5(a) and (b) shows that at higher thrusts the turbulent 
mixing effect is greater. Turbulent mixing reduces temperatures the most in regions 
where the fuel mass fraction is lowest. Therefore, the effect of turbulent mixing on the 
average fuel temperature is less than 1 percent. 

Figure 6 shows the effect of turbulent mixing on the average fuel temperature as a 
function of thrust for various pressures. Turbulent mixing affects the average fuel tem- 
perature by less than 1 percent for thrusts from 0. 445 to 4. 45 meganewtons (0. 1 to 
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Temperature, T, K 


6xl0 4 



8xl0 4 



Figure 5. - Radial temperature profiles at various axial positions for 
specific impulse of 1500 seconds and pressure of 5.05xl0 7 newtons 
per square meter (500 atm). Turbulence reduces temperature in 
downstream annular region. Temperatures near inlet and 
centerline are not very sensitive to turbulent effects. 
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.1 1 
Thrust, Mlb 

Figure 6. - Average fuel temperature as function of thrust for 
various pressure levels. Turbulence has very little effect 
on average fuel temperatures. Turbulence effect increases 
slightly with increasing pressure and is almost independent 
of thrust. 


1.0 Mlb) and for pressures from 10. 1 to 101.0 meganewtons per square meter (100 to 
1000 atm). The magnitude of the effect seems to be independent of thrust but slightly de- 
pendent on pressure. The authors conclude that in the range of conditions considered, 
heat transfer by turbulent mixing can be neglected in estimating the average fuel temper- 
ature. 

The average fuel temperatures in this report are about 15 percent higher than those 
given in reference 8. That reference considered a fueled region with a larger length-to- 
diameter ratio; temperatures were calculated only at the midplane of the fueled region. 
The axially averaged fuel temperatures in this report correspond more closely to the 
step fuel distribution, rather than the Gaussian fuel distribution of reference 8. This 
correspondence is partly due to the fact that the consideration in this report of buffer 
layer, wall coolant, and bypass flow resulted in reduced mixing. 


CONCLUSIONS 

From this study of heat transfer in a coaxial-flow, gas-fueled nuclear rocket with a 

specific impulse of 1500 seconds and a thrust between 0. 445 and 4. 45 meganewtons (0. 1 

7 

and 1. 0 Mlb) at an operating pressure of 5. 05x10 newtons per square meter (500 atm) 
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and for a dimensionless turbulent viscosity e + of 1000, the following conclusions were 
drawn: 

1. Radiation is the dominant mode of heat transfer for temperatures above 25 000 K, 
and turbulent mixing is the dominant mode below 20 000 K. 

2. For temperatures greater than 50 000 K, effective radiation conductivities are 
several orders of magnitude greater than effective turbulent conductivities. 

3. The heat-transfer process is not a constant-properties problem. For typical re- 
actor heat fluxes and temperature gradients the necessary effective conductivity corre- 
sponds to fuel temperatures of about 63 000 K and propellant temperatures of about 

35 000 K. 

4. The effective conductivity of a mixture of fuel and propellant is a strong function 
of the properties of the pure propellant and a lesser function of the properties of the pure 
fuel. 

For variable eddy viscosity, a function of both radial and axial position, the effects 
of turbulent mixing on the local and average temperatures are as follows: 

5. Near the inlet and the centerline of the reactor, the local temperature is reduced 
by less than 1 percent when turbulent mixing is included. 

6. The maximum change in local temperature at a thrust of 0.445 meganewton 
(0. 1 Mlb) was approximately 30 percent, while at 10 times this thrust (4.45 MN 
(1. 0 Mlb)) the change was equal to about 90 percent. 

7. The change in average fuel temperature is about 1 percent for all thrust levels. 

8. The change in average fuel temperatures is not very sensitive to pressure levels. 
It ranges from about 1/4 percent at 1. 07xl0 7 newtons per square meter (100 atm), while 
at 10 times this pressure it was about 1 percent. 

The preceding conclusions were drawn for a particular reactor geometry and flow 
distribution. Although these choices could affect eddy viscosities, the changes would not 
be by orders of magnitude. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, August 6, 1968, 

122-28-02-33-22. 
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APPENDIX A 


NUMERICAL ANALYSIS 


This appendix describes the numerical methods used to solve the equations derived 
in the heat -transfer analysis section. These equations are the energy equation and the 
equations defining average temperature, average density, average source strength, and 
the variation of the source strength. In these equations the fuel-propellant boundary 
temperature and heat flux are known boundary conditions. 

The energy equation is a nonlinear integral equation and must be solved by an iter- 
ative process. In the iterative process a temperature distribution was assumed. Based 
on this assumption the average quantities of temperature, density, and source strength 
were calculated. From these average quantities a new temperature distribution was cal- 
culated. This process was repeated until the average quantities changed by less than 
0.01 percent. 

The initial temperature -distribution calculation was based on an initially uniform 
source -strength distribution. For this source -strength distribution (p^/p^ s 1), the 
average source strength is 


Q 



(Al) 


For this average source strength the heat flux is 


q = 



(A2) 


When the heat flux and Fourier's conduction law are known, a Range -Kutta method can be 
used to calculate an initial guess for temperature distribution from the outer edge inward 
to the center line. If a uniform grid is set up in the radial direction r with Ar being 
the grid spacing, an approximate relation can be set up between adjacent grid points: 


<*n = 


-K n (T 


n+1 

Ar 


- T n> 


(A3) 


Equation (A3) can be solved by using a modified Newton-Raphson technique. For numeri- 
cal purposes equation (A3) was rewritten in the following forms: 


0 = In 


-Wl - T n> 


Arq. 


n 


(A4) 
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(A5) 


w = J Tfl - V 1 . 

V T n+1 

Equations (A4) and (A5) were used because of large property variations. For this tem- 
perature distribution the average temperature was calculated by using the following defi- 
nition: 




P f (T,y)r 


dr 


(A6) 


In order to calculate the average temperature, the Newton-Raphson technique was again 
used. 

For this average temperature the average density was calculated as 

p f = p f ( T, l) (A7) 


With this average density the average source strength was calculated as 


Q = 


r e^e^f 



(A8) 


With the calculated distribution and the calculated average quantities used as an initial 
guess, a Newton-Raphson technique was used to solve for the correct value of the aver- 
age source strength. The following equation was used in the Newton-Raphson technique 



- Q 


T TT 

guess 


(A9) 


where Qgyggg is the Q ca j from the previous iteration and Q ca j is calculated by using 
equation (A10) instead of (A2). 


q(r) = 




(A10) 


This process is repeated using equations (A10) and (A4) to (A8) until a convergence cri- 
terion is met. The convergence criterion used in all the modified Newton-Raphson tech- 
niques was based on the absolute value of the relative change in the independent variable 
being less than some specified quantity. 
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The mass fraction was obtained from reference 18 by using a linear interpolation 
method. The mass fraction is a function of the radial position r. The value for eddy 
diffusivity is also obtained from reference 18 in the same manner. 

The properties used to calculate conductivity and the density of the fuel are computed 
from tables. The properties for both fuel and propellant are viscosity, specific heat, 
Rosseland mean absorption coefficient, and molecular weight. These properties are a 
function of temperature and were computed by log-log interpolation. 
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APPENDIX B 


SYMBOLS 


Rosseland mean absorption co- 

Ar 

increment in radial position, m 

efficient, 1/m 

r e 

fuel-propellant boundary, 

specific heat, j/(kg)(K) 


y = 0.01 

thrust, N 

T 

temperature, K 

gravitational acceleration, 

U 

axial velocity, m/ sec 

9. 8 m/sec 

W 

function defined in eq. (A5) 

enthalpy, J/kg 

X 

mole fraction of fuel 

containment factor, defined in 
eq. (15) 

y 

mass fraction of fuel 

specific impulse (thrust divided 

z 

axial position, m 

e 

eddy diffusivity of momentum, 

by weight flow rate where ac- 


m 2 /sec 

celeration of gravity is 


o 

9. 8 m/sec ), sec 

e + 

ratio of eddy diffusivity to mo- 

2 

conductivity, J-m/(m )(sec)(K) 


lecular kinematic viscosity 

length, m 

V 

ratio of actual temperature to 


blackbody temperature at fuel 

molecular weight 


propellant boundary 

mass flow rate, kg/sec 


ratio of fuel flow rate to total 

3 

particle density, 1/m 


flow rate 

2 

pressure, N/m 
power, J/sec 

^noz 

ratio of actual I g p to ideal I g p 
for a nozzle 

2 

volumetric heat generation, 

P 

viscosity, (N)(sec)/m 

J/(m 3 )(sec) 

V 

2 

kinematic viscosity, m /sec 

heat flux, J/(m 2 )(sec) 

P 

O 

mass density, kg/m 

gas constant, 8. 3143 J/(K)(mole) 
Reynolds number 

a 

Stefan-Boltzmann constant, 
5. 67X10 -8 W/(m 2 )(K 4 ) 

radius, m 

radial dummy variable of inte- 

<P 

function defined in eqs. (8) 
and (9) 


gration 



I 


I I 


i mil ii ii i mini i iiiiiiii min in in ■iiiiiiii 


inn i hi im mi n ■ hi ami mi 


Subscripts: 

n 

dummy index referring to radial 

b 

buffer layer 


position 

cal 

calculated on present iteration 

P 

propellant 

e 

edge of fuel 

rad 

radiation 

ex 

exhaust 

t 

total 

f 

fuel 

turb 

turbulence 

guess 

guess of calculation on present 

w 

cavity wall 


iteration 

Oi 

dummy index referring to fuel 

h 

heat 


or propellant 

j 

jet, fuel stream at inlet condi- 

Superscripts: 


tions 


average 

i 

dummy index referring to spe- 

* 

blackbody conditions 


cies summed over 



m 

momentum 



mod 

moderator 
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TABLE I. - OPERATING CONDITIONS 


Parameter 

Range 

Value used 

Total reactor power, MW 

3970 to 39 700 

19 850 

Specific impulse, I g p, sec 


1500 

Thrust, F, MN 

0. 445 to 4. 45 

2.22 

Fuel- to total -flow -rate ratio, 7 ^ 


1/31 

Total mass flow rate, m ex , kg/ sec 

34 to 340 

170 

Fuel mass flow rate, m f , kg/ sec 

1. 1 to 11.0 

5.5 

Propellant mass flow rate, m^, kg/sec 

32.9 to 329 

164 

Pressure, P, MN/m^ 

10. 1 to 101 

50.5 

Nozzle efficiency, 7? noz 


0.90 

Propellant inlet temperature, K 


1600 

Exhaust enthalpy, h ex , MJ/kg 


133 

Ratio of propellant to buffer layer 

— 

5 

velocity, U p /U b 



Ratio of propellant bypassed to total 


2/3 

propellant flow 



Ratio of propellant used for wall cooling 


1/6 

to total propellant flow 



Ratio of injected propellant velocity to 


25 

fuel velocity 



Cavity diameter, m 


3 

Cavity length, m 


2. 1 

Inlet of fuel region, m 


2. 1 

Inlet of buffer layer, m 


2.55 

Total heat dissipated in moderator. 


10 

percent 
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